Phase Behavior in Thin Films of Confined Colloid-Polymer Mixtures 
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Using self-consistent field and density-functional theories, we first investigate colloidal self- 
assembling of colloid/polymer films confined between two soft surfaces grafted by polymers. With 
the increase of colloidal concentrations, the film undergoes a series of transitions from disordered 
liquid — > sparse square — + hexagonal (or mixed square-hexagonal) — > dense square — -> cylindric 
structures in plane, which results from the competition between the entropic elasticity of polymer 
brushes and the steric packing effect of colloidal particles. A phase diagram displays the stable 
regions of different in-layer ordering structures as the colloidal concentration is varied and layering 
transitions as the polymer-grafted density is decreased. Our results show a new control mechanism 
to stabilize the ordering of structures within the films. 
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Introduction 

Dispersions of colloidal particles in polymer solutions 
form the basic ingredients of a wide variety of systems 
which are essential to life and major industries) 1 ^ In 
particular, assembling colloidal nanoparticles into or- 
dered structureo 3 ! 4 ^ 6 ^ 8 ! 9 . 10 ! 11 . 12 ! 13 ! 14 ! 15 ! 16 ! 17 ! 18 is impor- 
tant to the fabrication of nanomaterials with high per- 
formance in size-dependent optical, electrical, and me- 
chanical properties. A major challenge in this field is 
to assemble monodispersed colloids into highly ordered 
structures with well-controlled sizes and shapes. Usu- 
ally, the dispersed structures depend upon the effective 
interaction between the colloidal particles themselves and 
between the colloidal particles and their surroundings. 
For example, the interaction between colloidal particles 
in a solution of nonadsorbing polymers is influenced by 
the entropic excluded-volume (depletion) forces which 
were first recognized theoretically by Askura and Oo- 
sawa (AO)fi 9 - and thus a monodispersed colloidal system 
may typically crystallize into a close-packed structure or 
body-centered cubic lattice with high symmetry. 

Thin films made of mixtures of polymer and colloidal 
particles can exhibit interesting structures, which are 
particularly useful for nano (bio) technology ranging from 
nanopatterning/lithography to changes of surface prop- 
erties. Surface confinement modifies the bulk behavior of 
the system by breaking the symmetry of the structure, 
and therefore can be used to make low-dimensional order- 
ing films. For nano-sized colloidal particles, the confine- 
ment of the suspension is a straightforward approach for 
the formation of colloidal crystals. For example, when 
the colloid/polymer mixtures are sandwiched between 
two parallel hard surfaces, layered structures may form, 
but the in-plane structure of the particles is always uni- 
form or close-packed. Polymer brushes grafted to sur- 
faces can modify the properties of the surface^S* 21 * 2 ^ 
which has many applications in colloidal stabilization 
and bio-compatibility. Interestingly, the grafted poly- 
mer chains may lead to a long-ranged repulsive force be- 
tween the particles, which stabilizes the ordering of the 



nano-sized particles and controls the symmetry of the 
structure formed in a desired direction. Such a control is 
required for a wide range of applications using structure- 
dependent optical, electrical, and mechanical properties 
of the colloidal crystals. It has been noted 23 i 24 i 25 i 26 i 27 
that thin films made of charged colloidal suspensions have 
various layered crystalline states if the film is confined 
between two parallel plates which repel the colloidal par- 
ticles. Here, we show that for colloid/polymer mixtures, 
using polymer brush surfaces is an effective way of chang- 
ing the in-plane structure of the colloidal film. 

The purpose of the present article is to demonstrate 
a novel route of spatially organizing the colloid disper- 
sions in polymer solutions confined between two polymer- 
grafted plates. We find that with the deformation of 
grafted chains by the particles, the deformed brushes 
can exert an additional interaction between nano-sized 
colloidal particles, which stabilizes the formation of low 
symmetric ordering structures of colloids. Such a brush- 
mediated colloidal effective interaction will compete with 
steric (center-force) effects of colloids, which may gen- 
erate structures with more complicated lattice symme- 
tries and undergo the symmetry-changing transitions in 
plane. In the present model, layering structures are 
formed due to confinement, and more importantly, the 
in-layer structural transitions are observed, depending 
upon the colloidal concentration and the polymer-grafted 
density. The results are summarized in a phase diagram 
displaying the stability regions of self-assembled order- 
ing phases as functions of the colloidal concentration and 
the polymer-grafted density. The proposed approach is 
physically achievable for the generation of novel ordered 
structures. A successful example indicating the impor- 
tance of such a self-assembling mechanism is the recent 
experimental verificationiS of the predicated phase be- 
havior of the copolymer/particle mixtures in confined 
geometries, 28,29 Another motivation of the present work 
is due to the fact that although many experiments have 
demonstrated the success of fabricating self-assembly 
nano-structures, the understanding of the physical mech- 
anism that drives the self-assembly of the colloidal par- 
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tides in confined geometries is still lackingiS^SLiffliSI 



Model Details 

We consider a colloid suspension confined between two 
planar plates separated by a distance L z along the z- 
axis. The two plates which are grafted with ribr A- type 
homopolymer chains, are horizontally placed in icy-plane, 
and mechanically fixed at z = and z = L z , respectively. 
All B-type free homopolymer chains in solutions are flexi- 
ble with the same polymerization N and statistical length 
a, and incompressible with a segment volume p$ 1 as in 
A- type grafted chains. The volume of the system V is 
L z , where L x and L y are the lateral lengthes 
of surfaces. The grafting density is defined as a = rib/ (2 X 
L x x L y ); the average volume fraction of grafted chains 
is (f)b = n.bN Pq 1 /V , the colloid (j> c = (1 — 4>b)tpc where V'c 
is the colloidal concentration of the colloid-polymer mix- 
ture, and the free polymer 4> p — 1 — 4>b — 4>c- Recently, 
the self-consistent field theory (SCFT) has been proven 
to be powerful in calculating equilibrium morphologies 
in polymeric systems, 32 i 33 i 34 i 35 i 36 i 37 i 38 i 39 i 40 i 41 i 42 while col- 
loidal particles can be treated by density-functional the- 
ory (DFX)iii 4 ^i 4 ^ 4 i to account for steric packing ef- 
fects of particles. In our calculation, we use the hybrid 
SCFT/DFT approach developed in refs 41 and 42, and 
the free energy F for the present system is given by 



Wi{r) qi , and = -^V\+ + W<(r)g+ . For the 



NF 



p k B TV 



-4>b ln( 



■)-<Mn(^ 



Vct> p 



v 



dr[xbpNcp b (r)<p p (r) + XbcNcp b (r)(p c (r) 



+X P cN(pp(r)tp c (r) - Wb{r)ifb(r) - W p {v)ip p (r) 
-W c {v) Pc {r) - e(r)(l - Mr) ~ <p P (r) ~ ^W) 
+p c ^ hs (7p c )} , (I) 

where fcs is the Boltzmann constant, and T is the tem- 
perature, a is the volume ratio of the colloidal particle of 
radius R to polymer chain: a = VR ^ a , where vr — ^irR 3 . 
y&(r), <p P {r), and <p c {r) are the local volume fractions of 
brushes, free polymers, and colloids, respectively. £(r) 
is the potential field that ensures the incompressibility 
of the system for dense colloid-polymer mixtures. p c 
stands for the particle center distribution, and the lo- 
cal particle volume fraction is then given by </9 c (r) — 



grafted chains in the field W b (r), the initial condition is 
q\(x,y,z = or L z ,0) = 1, q\(x,y,z ^ or L z ,0) = 0, 
and qi(x, y, z, 1) = 1, which means that the end of brush 
chains can move on the plates, although the total num- 
ber of chains on surfaces is fixed. These have been called 
liquid brushes, in contrast to solid brushes where the im- 
mobile chains are anchored onto the surfaces^ Here we 
consider the liquid brush case because of its wide-ranging 
applications in colloidal and biological systems. For free 
polymers, the initial condition is 92(2;, y, z, 0) = 1. The 
last term in Eq.(l) is DFT termi 3 . accounting for the 
steric interaction between particles, and the excess free 
energy ^(5? c ) per particle is from the Carnahan-Starling 
function 4 ^ ^hs(x) = 4 ^~ 3 J 2 , with the weighted particle 
density Tp c (r), 41 Tp c {r) = ^ J {r , l<2R dr'p c (r + r'), where 
V2R is the volume of a sphere of radius 2R. 

By minimizing the free energy in Eq. (I) with respect 
to W b {r), W p (r), W c (r), <p b (r), ip p (r), ip c {r), and £(r), 
we obtain a set of self-consistent equations describing 
the equilibrium morphology of films, which can numer- 
ically be solved by the real space combinatorial screen- 
ing algorithm of Drolet and Fredricksoni22ii& The free 
energy minimization of the system with respect to the 
selected simulation sizes is performed. 46 We fix N = 100, 
L z = 50a, and the colloid diameter 2R = 0.6i?o, where 
i?o = aN? characterizes the natural size of free poly- 
mers (typically 10 ~ 100 nm). On the other hand, based 
on the assumption that the interaction between differ- 
ent chemical species should be small enough and poly- 
mer brushes are chemically neutral to free polymers and 
colloids, we choose the Flory-Huggins interaction param- 
eters XpcN = 5.0 and XbpN — XbcN — 8.0. All the sizes 
are in units of a. 



vr Jir'|<-R dv' p c (r + r')i 41 i 42 The x's are the Flory inter- 
action parameters between the different chemical species. 
Q b = J dvqi (r,s) qf (r,s) and Q p = j drq 2 (r,l) are 
single-chain partition functions for brushes and free poly- 
mers under the self-consistent fields Wb (r) and W p (r) , re- 
spectively. Q c = J dr exp[— W c (r)] is the partition func- 
tion for colloids under the self-consistent field W c (r). The 
end-segment distribution functions qi(r, s) and gf~(r,s) 
represent the probability of finding monomer s at posi- 
tion r respectively from two distinct ends of chains, which 
satisfy modified diffusion equations = 2 -^-V 2 qi — 



Results and Discussion 

We first explore the formation of ordered structures 
of colloidal crystallization in confined colloidal solution 
film. Figure 1 shows a series of particle-dispersed mor- 
phologies of varying the colloidal concentration ip c , and 
Figure 2 is the top views of profiles for polymer brushes 
grafted to one of the plates. We observe the forma- 
tion of two layers of colloidal particles which are close 
to the polymer brush. When the colloidal concentration 
is small (for example, ip c — 0.29, as in Figure 1(a)), col- 
loids distribute uniformly in plane, and the brush-formed 
interfaces remain level(Figure 2(a)). With increasing the 
colloidal concentration, colloidal crystallization may ap- 
pear, and inhomogeneous colloidal distribution will lead 
to the deformation of polymer brushes: caving where 
colloids exist and protruding where not. The deformed 
brushes prefer to stabilize the formation of low symmet- 
ric ordering structures of colloids, and with increasing 
the colloidal concentration, the particles grow anisotrop- 
ically along one-dimensional direction^ We see from 
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FIG. 1: Distributions of colloidal particles with increasing the colloidal concentration for the brush density a — 0.2 and 
L x — L y — 50. The colloidal particles between two layers are alternatively packed each other, (a) ip c = 0.29, (b) ip c = 0.30, (c) 
ip c = 0.35, (d) Vc = 0.38, and (e) Vc = 0.45. 




FIG. 2: Top views of the density profiles of polymer brushes. Parameters are the same as in Figure 1. 



Figures l(b)-l(e) that for different colloidal concentra- 
tions, colloids can assemble into sparse square, hexag- 
onal, dense square, and cylindric structures, depending 
on the competition between the entropic effect of poly- 
mer brushes and steric packing effect of the particles. 
Figures 2(b)-2(e) show the corresponding shapes of de- 
formed brush surface due to respective colloidal struc- 
tures. When i(; c = 0.30, the average distance between 
colloidal particles is relatively large, and thus the brush- 
mediated interaction between colloids dominates over the 
steric packing effect of colloids, leading to in-layer square 
structures for possibly large conformational entropy of 
polymer brushes (Figure 2(b)). When ip c = 0.35, the sta- 
ble state becomes hexagonal one(Figure 1(c) and 2(c)), 
indicating that colloidal steric packing effects dominate. 
However, a dense square structure appears at ip c = 0.38 
(Figure 1(d)), due to the requirement of entropic elastic 
release of deformed brushes. Further increase of colloids 
greatly limits the conformational space of brushes, and 
the entropy of brushes may be released by the formation 
of cylinder arrangement of colloids, because in this case, 
the free ends of grafted chains can at least move freely 
between parallel cylinders, as shown in Figure 1(e) and 
2(e) when ip c = 0.45. Figure 2 shows that the polymer 
brush is soft, and its density variation changes with the 
colloidal concentration. The range of the density vari- 
ation in the direction perpendicular to the substrate is 
comparable to the size of the colloidal particles. The 
brush density close to the substrate remains unchanged 
with varying colloidal concentrations for a given grafting 
density. 

Figure 3(a) shows the entropy of a grafted chain as a 
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FIG. 3: (a) A grafted-chain entropy (-it 11 ) vs the colloidal con- 
centration; (b) the averaged steric energy per colloidal particle 
( ) vs the colloidal concentration. 



function of the colloidal concentration ip c . We set the 
grafted chain's entropy |kt = 0, when the dispersion 
of colloidal particles is uniform in plane {tp c = 0.29). 
The entropy of polymer brushes, with colloidal concen- 
trations varied from ip c — 0.29 to ip c — 0.32, increases 
with the formation of colloidal square structures. How- 
ever, the brush entropy curve suddenly decreases with the 
influence of steric interaction of colloids, and a hexag- 
onal structure is formed when tp c = 0.35. With fur- 
ther increase of colloidal concentrations, the deformed 
brushes begin to dominate over the steric energy of col- 
loids again, which prefers to release the conformational 
space of brushes with the gradual formation of square 
structures. The brush entropy arrives at the maximum 
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at ip c = 0.38, corresponding to the formation of dense 
square lattice structure (Figure 1(d)). After that, brush 
entropy decreases with the increase of colloidal concen- 
trations, and colloidal cylinder structures begin to form. 
The complete cylindric structure appears at ip c = 0.42, 
and remains unchanged even for relatively larger col- 
loidal concentrations. Usually, the entropic gain or loss 
due to the conformational change of polymer chains may 
signify the structural transitions of the polymeric sys- 
tem. For example, Balazs and co-workersM^ studied 
an AB diblock-nanoparticle composite, in which the par- 
ticles and diblock chains form spatially periodic struc- 
tures. They observed a similar behavior of A-block en- 
tropic free energy when the volume fraction of the parti- 
cles is increased^ In their work, however, the conforma- 
tional entropy of the A-block only shows the stretching 
and compressing of the A block chains along the chain 
direction, because they considered the formation of spa- 
tially periodic structures in two-dimensional copolymcr- 
nanoparticle composites in the plane of the polymer 
chains. On the contrary, in the present three-dimensional 
confined mixtures of colloidal particles and free homo- 
polymer chains, we emphasis a control mechanism in 
the lateral direction of the colloidal suspension driven 
by the deformed chains in the polymer brush. Actu- 
ally, the entropic free energy in Figure 3(a) depicts the 
shape changes of brush surface (Figure 2) during the vari- 
ation of the colloidal concentration for a given grafting 
density a = 0.2. The configurational entropy contribu- 
tion of the grafted chains originates from the free ends 
of the brushes comparable in size to the colloidal parti- 
cles, since the bottom part of the chains close to the sub- 
strate is strongly stretched and remains nearly unchanged 
with varying colloidal concentration for a fixed grafting 
density. This means that the entropic effect of the de- 
formed brushes by the particles comes mainly from the 
lateral conformation fluctuations of the free chain ends, 
and that the averaged conformational entropy contribu- 
tion of the polymer brushes along the chain-stretching 
direction can be ignored. Such an entropic effect com- 
petes with the steric packing effect of the particles in the 
plane, leading to structure changes in the plane of the 
colloidal particles. Here, a significant advantage of using 
soft polymer-grafted substrate is that the deformed brush 
provides an in-plane long-ranged anisotropic interaction 
between the colloidal particles. In Figure 3(b), the av- 
erage steric energy per colloidal particle varies with the 
increase of colloidal concentration, which clearly reflects 
the excluded- volume effect of colloids. As a whole, steric 
energy increases with increasing the colloidal concentra- 
tion ip c . Interestingly, we find that the slope of the curve 
has apparent difference in different ranges of ip c . Particu- 
larly in the region ip c ~ 0.32 — 0.38 where the morphology 
inclines to adopt hexagonal or mixed square-hexagonal 
structure, the slope is distinctly lower than that of other 
regions. This means that the steric energy can be re- 
leased by forming hexagonal or mixed hexagonal/square 
structures. 



Further, we also observe the appearance of layering 
transitions with the increase of effective film thickness 
via decreasing the height of polymer brushes. The ef- 
fective thickness of film can be defined by subtracting 
the brush height from the distance L z between two sub- 
strates, i.e., d e ff = L z — 2h, where the dry brush height 
h satisfies h — crN/po in the incompressible system of 
brush and homopolymer molecules4£ Thus, the effec- 
tive thickness d e /f of the film is changeable, and in- 
creases with decreasing the polymer-grafted density a. 
Figures 4(a)-4(c) show the possible three-layer structures 
such as mixed sphere-cylinder, complete cylinder, and 
even sphere structures for the selected values of a and 
ip c . To clarify the phase stability of different ordering 
structures^ we calculate the phase diagram as functions 
of the grafting density and the colloidal concentration, 
as shown in Figure 5. Here, we choose the grafting den- 
sity a within the range [0.15, 0.21] for two reasons: on 
the one hand, the entropic change of deformed brushes 
should be large enough to cause effective brush-mediated 
interaction between colloids, requiring that the value of 
<r can not be too small; on the other hand, the effective 
film thickness should be thick enough (namely a is rea- 
sonably small) to ensure the possibility of colloidal crys- 
tallization. We see from Figure 5 that when the graft- 
ing density is relatively large, the two-layered colloidal 
structure undergoes disordered liquid, square, and cylin- 
der ones. In this case, the brush-mediated interaction 
always dominates, which drives the self-assembly of col- 
loids into square lattice at small ip c and cylinder structure 
at larger ip c . The hexagonal phase disappears at the mid- 
dle concentrations of colloids. For intermediate values of 
polymer-grafted densities (a — 0.21 ~ 0.18), two-layered 
structure remains unchanged, but a reentrant structure 
transition between sparse square — > hexagonal (or square- 
hexagonal) — » dense square lattice structures appears as 
a result of the competition between brush-mediated in- 
teraction and steric packing effects of colloids. When 
a is less than 0.18, layering transition may appear for 
sufficiently large values of ip c . The three-layered sphere- 
cylinder structure (CS3 phase) first appears, where cylin- 
der structures closing to brush surfaces are formed be- 
cause of colloidal particles favoring the aggregation onto 
brush surfaces at first, and the middle layer still retains 
discrete spherical structure. With increasing the colloidal 
concentration, the cylinder- forming monolayer becomes a 
template which drives the ordering of other colloidal par- 
ticles in the middle of the film with the help of the defor- 
mation of free polymers in the suspension, and then the 
three-layered cylinder structure (C3) appears. If a is de- 
creased below 0.165, three-layered spherical structure (S3 
phase) may be formed for intermediate ranges of ip c , be- 
fore the appearance of mixed sphere-cylinder structures. 
The phase diagram shown in Fig. 5 also provides useful 
information about other structural changes of the system. 
For example, one can observe the structural and layering 
transitions with increasing effective film thickness (equiv- 
alently expanding the distance between the two plates) 



FIG. 4: The formation of three-layer structures: (a) sphere-cylinder mixture when a = 0.17 and ip c — 0.36; (b) cylinder 
structure when a = 0.17 and tp c = 0.42; (c) spherical structure when a = 0.16 and tp c = 0.32. 




FIG. 5: The phase diagram as a function of the grafting den- 
sity and the colloidal concentration. £2: two- layered sparse 
or dense square structure; M2: two-layered hexagonal or 
mixed square-hexagonal structure; C2: two- layered cylindric 
structure; S3: three-layered spherical structure; CS3: three- 
layered cylindric-spherical mixed structure; C3: three-layered 
cylindric structure. 



Conclusions 

We study the phase behavior in thin films of confined 
colloid/polymer mixtures. Depending on the competition 
between brush-mediated interaction and steric packing 
effects of colloids, the colloidal self-assembly can expe- 
rience a series of in-layer symmetry-changing transitions 
between square lattice, hexagonal, and cylinder struc- 
tures, which can be induced by changing the colloidal 
concentration. Furthermore, the layering transitions can 
be obtained by changing the film thickness of the sys- 
tem. We find that the colloidal particles lying on the 
polymer brush form a monolayer, in which the deformed 
chains produce an in-plane long-ranged anisotropic inter- 
action, which controls the ordering of the colloidal parti- 
cles. This layer of ordered particles then becomes a tem- 
plate on which the other colloidal particles in the middle 
of the film grow into a large ordered structure with the 
help of free polymers in the region. The results may pro- 
vide a helpful guide for fabricating the functionally useful 
microstructures in materials science of complex liquids. 
The suggested approach is also suitable to understand- 
ing colloidal self-assembling on the other soft walls such 
as the packing of biomacromolecules confined within the 
bio- membranes. 



via decreasing grafting density a for fixed colloidal con- 
centration ip c . If one changes the distance between the 
two plates while increasing the colloidal concentration, 
the colloidal structural transitions can be obtained along 
a tilted direction from the left to right in Figure 5. 
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To search the stable phases of the film, we first ascertain 

the symmetries of possible structure formation, and then 

calculate the free energy of the film. Furthermore, the total 

free energy minimization of the system with respect to the 

selected simulation sizes is required by adjusting the lateral 
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